NACA TM 1343 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



i 





NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 


THE EXCITATION OF UNSTABLE PERTURBATIONS 
IN A LAMINAR FRICTION LAYER* 
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With the aid of the method of small oscillations which was used 
successfully in the investigation of the stability of laminar velocity 
distributions in the presence of two-dimensional perturbations, the 
excitation of the unstable perturbations for the Hartree velocity 
distributions occurring in plane boundary -layer flow for decreasing and 
increasing pressure is calculated as a supplement to a former report. 
The results of this investigation are to make a contribution toward 
calculation of the transition point on cylindrical bodies. 
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I. STATEMENT OF THE PROBLEM 


If one wants to make a theoretical calculation of the profile drag 
of bodies in a flow for a certain direction of air flow, one must know - 
in addition to the pressure distribution - the position of the transition 
zone where the laminar boundary layer becomes turbulent. The separation 
point of the laminar layer forms a rearward limit for the transition 
point on a section of this body surrounded by the flow. It lies in the 
region of the pressure increase at the point where the velocity distri- 
bution in the boundary layer has the wall shearing stress zero. This 
separation point is a fixed point of the profile in the flow, the 
position of which does not shift due to a variation of the Reynolds 

number Re = (Uqo = velocity of air flow, t = chord of the body). 


One may calculate it according to the well-known approximation method 
of Pohlhausen which, for prescribed pressure distribution, provides for 
every profile point a velocity distribution of the boundary layer with 
a certain form parameter \ (X = -12 separation). As forward limit for 
the transition point, one may take the stability limit of the laminar 
layer with respect to small two-dimensional perturbations which were 
calculated in references 2 and 3 according to the method developed by 
W. Tollmien (ref. 1). According to this method, there exists for every 
form of velocity distribution in the boundary layer a so-called critical 


/U 6*\ / 

Reynolds number = 1-^ — 1 (Ua = local potential velocity, 

\ ^ / cr 

5* = local displacement thickness) below which all perturbations are 
damped; its value increases as the form of the velocity profile becomes 


fuller (with increasing X) . Where the Reynolds number Re* 


Ua5^ 

V 


formed with the actual potential velocity and the actual displace- 

ment thickness 5* exceeds this critical Reynolds number, there begins 
the instability of the boundary layer. In contrast to the laminar 
separation point, this stability limit is therefore not fixed on the 
profile for a normal pressure distribution but travels forward toward 
the stagnation point on bodies in flows with increasing Reynolds 
number Re (see fig. 1). 


The actual transition point which lies between these two limits - 
stability limit and separation point - is known to likewise shift forward 
with increasing Reynolds number Re. That it lies only a certain 
distance behind the stability limit (as shown by a comparison of experi- 
mental transition points and theoretical stability limits) is plausible, 
because the excitation of the unstable perturbations starts only at this 
limit point of the stability and must obviously have attained a certain 
degree before the instability further downstream leads to the breakdown 
of the laminar flow configuration. 
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For this reason it seemed necessary to calculate, or at least to 
estimate, for the velocity distributions in the laminar boundary layer 
in the entire instability range of the latter, the excitation of the 
unstable perturbations as well, with the objective of making an improved 
calculation of the transition point possible. 

For the velocity distribution on the flat plate in longitudinal 
flow (Blasius profile), H. Schlichting (ref. 4) has already determined 
the excitation quantity as a function of perturbation frequency and 
Reynolds number in a part of the instability range; and for the velocity 
distributions in the region of pressure increase, W. Tollmien (ref. 5) 
has explained the behavior of the excitation (in first asymptotic 
approximation for very large Reynolds numbers) in a very general manner, 
neglecting the effects of internal friction. 


II. THE GENERAL DIFFERENTIAL EQUATION 
DESCRIBING THE PERTURBATION 


Since the bases for the method of stability investigation have been 
discussed in detail in an earlier report (ref. 3)^ we can refer to the 
results attained there. 


Let U(x,y) and V(x,y) be the tangential and normal components 
of a plane steady boundary- layer flow; let x denote the length of the 
arc and y the normal to the profile contour. The stream function of 
the two-dimensional perturbation motion which we superpose on this 
basic flow is assumed to be 


t(x,y,t) = 9(x,y)e^°'^’^‘^^^ = cp(x,y)e^^^e^^°'’^"^^^^ 
c = Cj. + icj^ = Cj^a Pj. = Cj,a 


( 1 ) 


here t denotes the time, a the spatial circular frequency of the 
perturbation, the real part c^ of c its phase velocity, and the 

imaginary part c^ a measure for its excitation (ci> 0) or damping 
(ci < 0); besides Pi is the logarithmic increment of the excitation 
of the perturbation amplitude, and Pj. the circular frequency in time 
of the perturbation. 
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If we substitute the motion originating by superposition of the 
boundary-layer flow U, V, and the perturbation motion 


u 


V 


= ^ = (pieia(x-ct) 

Si ( • m ^ 




> 

ia(x-ct) 


C2) 


into Navier -Stokes ' differential equations^ we obtain - as was proved 
in reference 3 in detail - the differential equation describing the 
perturbation in the form 


(U - c)fcp’* - a^cp) - U* 'cp = ' - 2a2cp»' + a^) 

\ / oRe* ^ ' 


Re* = 


oRe* 

Uab* 


^ (3) 




Here the prime ( ’ ) denotes differentiation with respect to the wall 
distance y; the velocities are referred to the local velocity Ua at 

the boundary -layer limit, and the wall distance as well as the wave 

length ~ referred to the local displacement thickness 6*. 

In order to avoid misunderstandings it should again be emphasized 
that - in spite of the assumption that U, V, and cp be functions of 
the arc length x - only the form of the local velocity distribution is 
decisive for the stability investigation as one recognizes from 
equation (3). The immediate effect of the pressure gradient, however, 
is negligibly small as is also the influence of the x-dependence of the 
perturbation amplitude cp. 

The boundary conditions of the differential equation (3) result 
from the condition that the perturbation velocities u, v vanish at 
the wall and that the friction effect at the outer boundary -layer limit 
(U'* = 0) has disappeared. With their aid, the calculation of the 
excitation of the unstable perturbations may be reduced to an eigenvalue 
problem, the solution of which is discussed in section V. 

We deal first with establishing the particular solutions of the 
differential equation (3), limiting ourselves to small values of the 
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excitation quantity thus, the general solution of the differential 

equation describing the perturbation (3) may be represented in the form 

k ^ / ic- \ 

^ = II CvV = ^ Cv^k — o>vj Uo' = (4) 

There cp^ denotes the particular solutions for c^ = 0 which are 

obtained in the calculation of the limiting curve of the instability 
range (neutral curve) in the a, Re*-plane (ref. 3)^ and the o)^ 
signifies additional functions for > 0. 

We turn first to calculating the integrals cp^*, ^ 2 * 

additional functions 0 ) 2 . 

III. THE SOLUTIONS cp^*, cp^* OF THE DIFFERENTIAL 

EQUATION DESCRIBING THE FRICTIONLESS 
PERTURBATION FOR FINITE EXCITATION 


If oRe* is assumed to be very large, the differential equation (3) 
is simplified to the so-called frictionless -perturbation equation 


(U - c)(9' ' - a2(p) - U"cp = 0 


(5) 


This differential equation has a pole of the first order at the location 

U = c = Cj. + ici to which we coordinate the point y^* of the complex 

y-plane. In the neighborhood of this singularity, a fundamental system 
may be easily indicated by series development. In order to establish 

the connection with the case of the purely real c treated before 

(see ref. 3)> we first give the relation between the complex y^* and 
the wall distance y^ of the "critical" layer U = cr. 

From 


U(yc) = Cr 


( 6 ) 
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and 


U^Yc*) = Cj. + ici (7) 

there follows 

u(yc*) - U(yc) = ici = ^Yc* - Yc)Uo’ + • • • (8) 

and, with limitation to the terms linear in Cj^, therefore 

y. = 7c * ^ (9) 

We shall now indicate the construction of the solutions cp-j^*, 

of equation (5) for those special velocitY distributions U(y) bY which 
\re shall approximate the laminar boundarY profiles of Hartree (ref. 6) 
for the calculation of the excitation of the unstable perturbations 
in the same manner as we did before in reference 3 for the calcula- 
tion of their critical ReYnolds numbers. 


(a) Binomial VelocitY Distribution 
(Pressure Decrease) 

In the region of the pressure decrease, we used the approximation 
function 


U = 1 - (a - y)" n = 2 , 3 ,^, . • • (10) 


where a denotes the coordinate of the point of junction to the potential 
velocity. 


With the new variable 
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and with 


aj_ = a(a - yc) 


( 12 ) 


the perturbation equation (5) for indifferent perturbations then reads 


[1 . (1 . n)^ (r 


d^ 


- + n(n - 1)(1 - yi)’^~^cp = 0 (I 3 ) 


Using the relations (9) and (11), we now introduce the complex variable 

* y ■ yc* yc - yc* 

yi =: — = — r:^ = yi (1^) 


a - Vc 


a - 


Uo'(a - yc) 


With the abbreviation 


= 


=i 


Uo'(a - yc) n(cr - 1) 


(15) 


the differential equation describing the frictionless perturbation for 
nondisappearing excitation then is transformed into the form 


1 - (i - y^*) - (1 - 


' d^cp* 2 * ' 
,^1 


n(n - 1) 


(l - yi*)’"'^ + (n - 2 )ifi(l - y2*Y'^ 


= 0 (16) 


in which, according to our presupposition, only the terms linear in 
f]_ have been taken into consideration. 
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If one writes equation (l6) after multiplication by 


yi 


*2 


1 - (l - yi*)” + nif^ 1 - ^1 - 


and after division in the form 


- 0 

^ dy/2 


the first solution cp^^* is given by 



with 


0)1 


= £ «vyi'' 

v=o 


The series coefficients e^* are obtained from the 3v* 
to the recursion formula 


v(v 


v=l 




* T 

eo = 1 


The solution for all y^* is found to be 


* - - * r ^ 

Tl 


<32 = ^2e^ 


* *v 


CPi 


a - Yc 


In : 


(17) 

( 18 ) 

(19) 

according 

( 20 ) 

* ( 21 ) 
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In this equation, the logarithms are, however, ambiguous; one must cut 
the complex y^^^-plane in such a manner that 

3 * JT 

- I « < arg < | 

If we put 

III Vyi*'' = IZ ^yi'' + YZ. K»yi'' (22) 

V=2 V=2 V=1 

cpg* niay be represented for y^^ = real part of y-j^* >0 at first in the 
form 

00 00 

cp2* = 1 + YZ YZ ^'yi^ 

V=2 V=1 


with 






if-| \ 

yi + = ^P2 - — T “2 


ic. 


yi 


Ur 


(23) 


cD2 = 


a - yc 




V=1 




w 

2^1 In y^ y~ g»y/ 
v=o 


( 24 ) 
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In the transformation of equation (21) into equation (23) the 
relation 


®i* " ®i(^ ■ ^^1) 


(25) 


is put to usej it is easily obtained from the perturbation equation (5) 
if U(y) and cp-[_* for = 0, and for Cj^ = 0 in the neighborhood 

of the singularity U = c are developed into a Taylor series, since 


and 


2e, = (a -■ yc)H0- 


uo 


(26) 


2e±* = (a - yc)- 


Uq" + ici- 






Uo- . ic^:^ 


Uo’ 


(a - yc)^^" 


Ua' 


1 - 


iCi| 

(Uo" 

Uo*' 

(Uo' 


Uo' 


Uo" 


Uo" , 

= (a - yc)-^(l - ifl) 
Uo' 


(27) 


If we express in equation (27) the exact value 

Uo'"' 


Uo' Uo’ 
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again by the approximate value 
that due to the small value of 


term for the value of 


Uo’ 


f^^ this step is motivated by the fact 

one deals here only with a correction 


For y^ < 0, cp^* represented by the expression 


-Pa* -1*1: 


+ i arg 


lC-5 


CP 2 - — 0 ^ 


^ ^ Jt 

- ^ ir < arg y* < ^ 


(28) 


with 


(Do = 


- yc 


V=1 


h^yiV + — - In |yi| 


a - yc^yi 


2e^ In 


yi 


T1 vi^' 

v=0 




la - yc v=0 


(29) 


The term with i arg y^* in equation (28) was obtained by 

H. Schlichting (ref. 4) and W. Tollmien (ref. 5) by a discussion of the 
general perturbation equation (3) in the neighborhood of the singularity 
U = Cj. + icj^ in a similar manner as the "transition substitution" in 

the critical layer for pirrely real c by taking the friction effect 
into consideration. 


Since we are, in the calculation of the transition substitution, 
concerned with a representation of e^^ of maximum accuracy, we shall 

replace ej_ by the expression (26) of the Taylor development of the 
exact Hartree profile. 
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For the linear and the parabolic approximate distribution 
[n = 1, n = 2 in equations (10) and (l6)J H. Schlichting (ref, 4) 
has already given the solutions 92*’ calculation was con- 

tinued for n = 3 and n = 4. The coefficients 

'g^^ by are compiled numerically in table 1. For the convergence of 
the power series^ the reflections made in reference 3 are valid. 


(b) Sinusoidal Velocity Distribution 
(Pressure Increase) 

In the region of pressure increase^ the Hartree velocity distri- 
butions in the boundary layer (see ref. 3) were approximated by 
formulation in terms of a sine formula introduced by W. Tollmien 
(ref. 5). ■ 


U = Ug + (1 



(30) 


where s denotes the wall distance of the inflection point. 
With the new real variable 


^2 = o 


« y - Yr 
2 a - s 


(31) 


and with 


02 = 


-(a - s)a 


(32) 


the perturbation equation (5) for the indifferent perturbations then 
read 


~i / d^cp 


|^ln(y2 - Y2 s) + sin 72^/^ - a2^cpj + sin - y^^y = 0 


(33) 


MCA TM 13^+3 


13 


where 72 g was put equal to 


y 2 s = ^ 


( 3 M 


We now introduce the complex variable 


y/ = £ y.^ = y2 ^ £ yc_^ 
2 2 a-s 2 a - s 


= I2- 


ic-i 


Uq' 2(a - s) 


(35) 


With the abbreviation 


fo = - 


Uq* 2(a - s) (1 - Us)cos ygg 


(36) 


the frictionless differential equation, describing the perturbation 
for non-vanishing excitation, then is transformed into 


sin ^72* - y2s) + sin 723 - if2|cos ^72* - 72s) - ^os y2g| 

+ |sin^2* - y2s)-if2 cos (73* - yps^ 9 * = 0 


- a2^cp'^ 


If (37) is written like equation (16) in the form 


72*2 ♦ 9* E Bi*72*1 - 0 

dy2*2 i=i 


(37) 


(38) 


the solutions cpp*, 92* determined by analogy with the equations 

(18), (23), and (28). One has 


,, , ,1* = 72* n - - 

2(a - s) V=0 


”1 ' >^2 *'■ ’'S 2 (a - s) 


IC^ 


^ ■ uv 


(39) 
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with 




-I 


v=o 




(4o) 


For Y 2 > 0, the solution cp2 obtained 

00 TTCp,* 

q> 2 * = 1 + ^ 4>y*y2*'' + 2ei* — In y 2 ^ 

V=2 


- 1 ir < arg y * < 1 

2 2 2 , 


(41) 


If one puts 


bv*y2*'^ = YZ ^ + 1^2 Y ~ ^^2^ 

V=2 v=2 v=l 


(iv2) 


CP 2 for Y 2 ^ ^ represented in the form 


Cp 2 * - 1 + Y2 + ^^2 ■*■ ^ 

V=2 V=1 \ SI 

^ gvy2'||^ln J2 + 


if2 \ 


sin y 2 g cos y 2 g 


.2(a - s) 


1 C 

= % ^ (U- 


U, 


0 


7 ^2 


( h 3) 


with 


0)o = 


^ 2(a - s) 


n h,y2v . i iyi /j. , iiLw : 

V=1 2(a - s) ( y2 sin ygs cos y2s 


2e, In yo 5 


g..yo 


V 


(44) 
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In transforming (4l) into (43), one makes use of the relation 


®1* - ®1(^ 


sin Y2s cos y2s 


which follows from the two equations 


■ 


a - s V 
Jt Un 


(a - s) ^O" + ici 


TI ' ' ' 

^0 


Uq- + iCi- 


(a - s) Up” ^ ici/Uo” Uo’” 

« Uo’ ■ Uo'luo’ ■ Uo” 


(a - 3) Up” / ^ if 2 

JT Uq ’ i sin ygg cos y 2 j 


cpg* for y 2 < 0, is represented by the expression 


rp2* = 1 + ^ (a _ 3) ^1*1 1*' ^2* I ^ ^ ¥ 2 * 


cp2 0)2 

Up' 


- I Jt < arg < J )(48) 
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with 


a3.^ = 


2 2(a - s) 

2^1 ln|y2lZI ^^2 

' ' V=0 


f V/ - 

^ “V 2 a - 3 ly 2 sin y^s cos y2s ) 


+ 2iJTe3^ 


- "fPi 


2(a - s) sin y^g cos yg^ 


v=0 


Svy2 


(49) 


For the calculation of the transition substitution, we shall 
replace the term as in section III (a), by the accurate value 

(equation (46)) of the Taylor development of the exact Hartree profile. 

In a comparison of the relations (24) and (44) or (27) and (47) or 
( 29 ) and ( 49 ) > it is striking that, in the expressions for e^^* and 

for the sinusoidal boundary -layer profiles in case of pressure increase, 
the product sin J 23 ^28 appears in the denominator of several 

terms. The sign of this product is negative when the inflection point 
of the velocity distribution lies more closely to the wall than the 
critical layer (723 ^ critical layer thus lies in the part of 

the velocity profile showing concave curvature. The sign of the expres- 
sions divided by sin y^^ cos y^^ then is the same as the sign of the 

corresponding terms in the solutions for the velocity distributions in 
case of pressure decrease which, as is known, have concave curvature at 
every wall distance. If, in contrast, the critical layer is located 
between point of inflection and wall ^y 23 >0^, then, in the part of the 

velocity profile having convex curvature the signs discussed are reversed. 
In that case, when the critical layer shifts to the point of inflection 
itself ^y 2 s = 0\, the behavior of e^^* and 0^2 is regular since then 


e 


* 

1 


iCi 


2(1 - Us) 



( 50 ) 
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and - 



The series coefficients ®v*^ 


Sv^ 


hy for the 


sinusoidal basic velocity are given numerically in table 2. For the 
convergence of the series developments, the explanations of the earlier 
report (ref. 3) are again valid. 


IV. THE FRICTION SOLUTION cp^* FOR 
FINITE EXCITATION 


Besides playing a role in the critical layer U = c, the friction 
of the fluid is of importance also in the neighborhood of the wall where 
it occasions two more solutions 7*4* general differential 

equation describing the perturbation (3). 


If we introduce the variable 



with 


e* = 


oRe* 



+ 


IC-i 



(53) 


(54) 
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we obtain from equation (3) in "the limiting process € — »0 for large 
Reynolds numbers the differential equation 


dti*'' dti*2 


(55) 


The solutions of this differential equation are just as independent of 
the form of the velocity distribution U(y) as was the case for the 
solutions cpj, 94 for the excitation zero (ref. 1)} they read 

corresponding to reference 2: 



* 1/2 ( 1 ),( 2 ) 

' ®l/3 



(56) 


where signifies the Eankel function of the first and second 

kind, respectively. 

Since the Hankel function of the first kind increases for large 
wall distance beyond all limits, it cannot be contained in the general 
integral (4) so that we there may put = 0. 


V. STATEMENT AND SOLUTION OF 
THE EIGENVALUE PROBLEM 


After having found the particular solutions of the differential 
equation describing the perturbation (3)> we now state the eigenvalue 
problem, (which results from the boundary conditions of the perturbation 
equation) for investigation of the excitation of the unstable pertur- 
bations. 

At the wall, the tangential as well as the normal component of the 
perturbation velocity disappear; thus one has 


Ci*cpi*w + C2**P2*w + *-3**T*3*w = 0 

(57) 

CiVw' + CgVw’ + C3*cp3\' = 0 

(58) 
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At the point of juction y = a to the region of constant velocity, 

U" is zero and therewith cp*-" - a^cp* = 0, that is cp* = e“*^y, so that 
the third boundary condition reads 


+ ^ 2 *^ 2 * e. = 0 


with 


4*.. = cp, * + acp * , V = 1, 2, 

V a a a^ . 


A term *1063 appear because the particular friction solution 

cpj* has already been damped at the point of junction. 

In order that the three homogeneous equations (57)^ (58)^ and (59) 
may have a solution different from zero, the determinant must disappear 


> 

e- 

^2 w 


m * ' 

'Pi W 

* ’ 

^P2 w ■ 



$ * 

0 

1 a 

2 a 


determinant yields the 

'PS w 

^2%h\ 

-cpi- 

cp * ' 
^3 w 

^2%'^la 

- 


or with Schlichting' s abbreviations 


i = . Eholl = Fho») - F* 

^10* 'P3*v' 3o* 
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1 - ^l\^2\ 

yc^ - ^l%'^2*a 


E*(a,cr,ci) 


( 64 ) 


finally 


F* = E*(a,cr^ci) 


(65) 


Therein 


^0 


* 



( 66 ) 


The complex equation (65) is equivalent to two real equations in which 
the parameters a, Re*, Cp, and c^ are contained. If one limits 

oneself to the case of indifferent perturbations (ci = O), one obtains 
from them after elimination of Cj. a relation between the Reynolds 

number Re* and the wave length ^ of the perturbation. That is the 

equation of the neutral curve or the curve of constant zero excitation 
by which the stable and unstable perturbation states are separated and 
the lowest Reynolds number of which is the so-called critical Reynolds 
number Re*^pj when this critical Reynolds number is exceeded, the 
excitation of the unstable perturbations begins. 

These neutral curves were determined in a previous report (ref. 3 ) 
for the velocity profiles for decreasing and for increasing pressure 
calculated by Hartree. 


In the present investigation, we shall assume the excitation 
quantity Cj^ to be different from zero and investigate the curves of 
constant excitation enclosed by this neutral curve. 


For this purpose, we shall solve the complex equation (65) in a 
somewhat simpler manner than Schlichting (ref. 4 ), by determining not 

5 Ci S Ci 

the differential quotients — — and — 


at the location of the 


S a S Re* 

neutral curve, but the curves of constant excitation Ci > 0 directly 
by a method similar to that applied for the neutral curve Cj^ = 


however, without making use of this curve itself. We consider first 
the left side of equation (65). 
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Since we intend to limit ourselves to small values of excitation, 
we develop F* in the form 

F* . F(no*) = F(rio) t , (67) 


where tIq as before in reference 1 is defined by the relation 



According to equations (9)^ (5^)> ( 66 ), and ( 68 ) one now has 



Thus equation ( 67 ) becomes 



” (to) 
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One sees therefore that the form of velocity distribution which does 
not enter into the exact solution according to equations (63) 

and (56) does appear explicitly in the development of with 

respect to Cj^. 

The differential quotient 

dF dFj. ^dFi 

driQ dT)o <iT]o 


was determined by graphic differentiation of the function F^T^Qj. The 

numerical values of its real and imaginary part are given numerically 
in table 3 . Since real and imaginary parts nowhere disappear simultane- 
ously, the function f(t)q*) will be free from singularities in region 

around the function F^tiq) and the development (67) will be thus per- 
missible. 


We now consider the right side of equation (65) which is defined 
by the relation ( 64 ). 

Remembering the splitting of the frictionless solutions 
cpy'*l[v = 1,2) in the form of equation ( 4 ), we first put 


with 


furthermore 


A * 

V a 

= 0 

va 

■ ^ 


= %a' 


^ * 

> 

II 

ici 

'Pv W 

- — ^w 
^0 

cp * ’ 
w 

'^vw 

ici 

" “~T 



Uq 


(71) 
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For the binomial velocity distribution U = 1 - (a - y)*^ the 

iS of the additional functions (o,, 

Va 

^v' ^ given by the equations 


values of the additional functions o)-, 03 ., and their derivatives 

Va y^ff 


^a = ZI Sv 

V=0 




2a 


a - yr 




V + 

v=l a - yc 


OJ- 


la = 




a - yc v=i 


VSv 


o>, 


2a 


1 2 ei / 2 cpn , 

vhv + ^ + cpia + o)ial 


(a - yc)^ v=l (a - yc)2\ a - y^ 


O) 


Iw ^^Iw 




a - yc 




2e^ In 


. 1 , 


0 ) 


Iw 


CDo * 
2w 


00 






— - 03 

Iw 1 

^ - yc 

/ 



2 ei 

'Plw / 

(a - yc)^ 

a - yc\ 


>(72) 


1 1 


(a - yc)' 

'Piw'(^^ - * is * - yc) 


yiw/ 


< t / _ > 


uo U - yc 
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For the sinusoidal velocity distribution 

u = Us + (1 - Us) 81 h/ 5 L:l^ J 

- s 2 

in contrast, the corresponding formulas read 




V , "e^cpia/ 1 In y2a 


2(a - b)[^ a - s \^y2a y2s <=08 y2g 

2e^ In y2a‘“in 


t JT 
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— ■■ ■ ; -— + In ygaO^ia 
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g-/ 2 w' 


“iw = 2 _ g-/ 2 ^ 

v=0 
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« I V V 1 y2w| \ 
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A 


If the critical point shifts to the inflection point of the velocity 
distribution = s, y2s = "the equations for ^2a' > 

^2w' (equation (7^)) are simplified to 


0 ), 


'2a 


2(a - s) 


V=1 


V2s 


"7>la 
2(a - s) 


In yga 


0 ), 


'2a 


“2v 


ff2 

4(a - s)2 


vh^yg^v-l + 


v=i 


^^la 

2(a - s)y2j 


+ cpia' In yga 


la 


2(a - s) 
«2 


TL V2w'' + ' In |y2w| 

V=1 2(a - s) ' ' 


ijr3, 


9lw 


> (1^) 


0), 


'2w 


i|(a - s)2 

ift3cpi^' 
U(a - s)^ 


U(a - s)2 

5 " 5irrfc- " -fiv' l^swi 


We now split *P2w> ^2w > into real and imaginary 

parts according to the formulas 


'^2w ■ -^1 ^®1 + iNj 


CP2^’ = Ag + iBg 0)2^' = Mg + iNg 


(75) 


then one obtains 


^ ”1 * '('i ■ ^ ”i) 




92 * w ’ - ^2 { 7 / ^2 ■*■ ^ 2 ^ 

uq \ '^0 y 


J 


( 76 ) 
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Then we put the expression (6i+) into the form 


E*(a,Ci.,Ci) = E(a,cr)[^ + Ci(zi + iz2^ 


where E, as previously defined in reference 1, is defined by 

1 92w<I'la - '?lw^2a 
E ■, ; 

yc T>2 w ^la ■ '^^Iw '^2a 


First, limiting oneself to .the terms linear in one has 

"P2*w‘^l''a ^l^la ^ ®l^ia) 


^ ®l®la - ^(^l^la ^l*la 


ic ^ 

Pl*w®2*a ~ ^lw^2a ‘ “ (‘“lw®2a + 


'^2*w' Va " Via i^(Vla + Via) + 


i B2$ia - -^(Vla + 


~ 'Plw'<I’2a - .^Kw’-^2a + 
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Hence, there results after a short intermediate calculation 


E* ~ E 


IC-i 


Uo yc 


A 

Uo’ 


(Ni<I»ia + Biflia) + i(«>lw^2a + 9lwfi2a ~ ^l^la ~ Mi<I>ia) _ 


^i^la - 'Plw'^2a + iUl'^la 


(U2'^la + B2%a) + i(*^lv'^a + q>lw’%a - ^2^1a ~ M2^1a) 


^2^1a ■ ^lw'^2a + i®2®la 


(80) 


We put for abbreviation 

^1‘^la + ®l%a = % 

^lw^2a + 'Plw^a ~ Ai^ia - Mi<I>ia = ni 
N2<I>ia + B2ftia = “2 

‘^lw'^2a ■*■ "^^Iv'^a " ^2^1a “ M2^1a = «2 


} (81) 


^I'^la - 'Plw’^2a = Kl 
^2^1a ■ 'Plw''^2a = ^ 
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Then there follows from equation (80) 


E ~ El - 


ici 


™1 ^’^l 

Uo'yc 

H — - 

Uo’ 

Ki + 

k. 


m2 + m2 


= E 1 - 


icj ^ °^lKi + j^iBj^ia _ 1112K2 + n2B2*^ia 
Uo'yc Uq' Ki^ + K2^ + 


iC;L i^l^i - nilBi^>ia H2K2 - 
Uo'|ki^ + ^ 2^ + B2^*^'la^ 

If we finally introduce the designations 

2 ^ _3^ miKi + n3_B]_<Ii3_a ni2K2 + n2B2<I’ia 

^ ■ Uq' Kj _2 + Bj_ 24 i ^2 ■ 


1 hiiKi - m-|Bit-|„ UpKp - 


Uo'yc U(-)' Ki^ + Bi 2$^^2 k 2^ + Bp^la^ 


the representation (77) is attained. 

By means of the equations (70), (77), and (83) we thus have divided 
the two functions F* and E* (a,Cr,ci) appearing in equation (65) into 
real and imaginary parts and can now graphically solve this equation. 

For this purpose, we plot for constant c^^ for several values y^ the 
imaginary part of the function F* 


i(f*J = I 


1 + Uq"Wi 


Uo’yc 3 Uo'/dqo 


against its real part 


k(f») . k|f(bo)] 


- 'lo .7-^1 * 
Uo ycv 


Uo 


Uo' 
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and plot in this polar diagram to the same values and y^,, the 

imaginary part of the function E*(a,Cr,ci) 

= I(E)(1 + CiZi) + CiZ2R(E) (86) 

against its real part 

r(e*) = R(E)(1 + CiZ^) - ciZ2l(E) (8?) 

The points of intersection of the functions F* and E* connected one 
to the other hy the same pairs of values Cr, Ci yield first the values 
t)q pertaining to the corresponding a-values, and then with (68) the 
cxirve a(Re*) of the constant excitation Cj^. This plotting and cal- 
culation then are repeated for other values Cj^. 


VI. RESULTS OF THE CALCULATION 


As immediate result of the rather extensive calculations, the polar 
diagrams with the curves E*(a,ci,cr) and F* = for several 

Hartree velocity distributions U(y) in the laminar boundary layer are 
represented in figures 2 to 21. These boundary -layer profiles belong 
to the special class of pressure distributions 

Ua = const, x® (88) 

and are rigorous solutions of Prandtl's boundary- layer equations which 
were obtained with the aid of a Bush apparatus (ref. 6). We use, 
according to Hartree, for characterization of these boundary -layer pro- 
files the parameter 

3 = (89) 

1 + m 

The profiles with m > 0, 3 > 0 occur in case of decreasing pressurej 
the profiles with m < 0, 3 < 0 in case of increasing pressure. The 
profile 3=1 is the exact Hiemenz profile at the stagnation point of 
a cylindrical body, the profile 3=0 is the exact Blasius profile at 
the flat plate in longitudinal flow, and the profile 3 = -O.I98 is a 
separation profile (wall shearing stress zero) . 


># 


We have investigated the excitation for the Hartree profiles p = 1, 0,6, 0.2, 0, -O.IO, 
-O .198 (the last, however, only for larger values of the Reynolds number, because of the uncer- 
tainty in case of small Re*); as in an earlier report (ref. 3)^ replaced them by binomial 
approximation distributions in case of decreasing pressure and by sinusoidal approximation 
distributions in case of increasing pressure, with the aid of the formulas 
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is the dimensionless wall distance used by Hartree and the coordinate 
a of the point of junction to the potential velocity Ug is connected 
with the displacement thickness by the relation 


— = - 3 

& kg 


(92) 


the quantity kg*(3) may be found plotted graphically in figure 2 of 
reference 3* 

The wall distance of the critical layer indicated in figures 2 

to 21 is obtained if in equation ( 91 ) y is replaced by yc- 

The following basic remarks should be made concerning the graphic 
solution of the eigenvalue problem in figures 2 to 21: 

Since Cj^ has been presupposed so small that only the variations 
of the particular solutions cpy* linear in Cj^ need to be taken into 
consideration compared to the solutions cpy for c^j^ = 0 , the curves E* 
and F* appear for equal c^-interval as "equidistant" curve families 

in the sense of equations (77) and (70). Actually, this "equidistance" 
will be lost in case of higher values of Cj^j however, the calculation 

expenditure would increase intolerably even if only, for instance, the 
terms quadratic in Ci were to be taken into consideration. In this 
sense, the curves 

3j^6* Cj^ab* 

”uT " Ua 


against a5* represented in figures 22 to 26 (from now on we use dimen- 
sional quantities) which result from the evaluation of the polar diagrams 
also are to be interpreted as approximations (in this and the following 
figures, the value which is, according to equation (l), physically 

more important and which characterizes the logarithmic increment of the 
excited perturbation amplitude has been plotted instead of c^). Here, 
as in Schlichting's report (ref. i^) mentioned before, only the derivative 


d(a6*)\Ua ) 
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at the location 


= 0, that is, the slope of the 


curve at 


the a6*-axis is rigorously correct. H. Schlichting interpolated 
between the base points at which he had determined the slope directly, 
although with a much higher calculation expenditure, with a curve of the 
third degree 


= aQ + aj^(a&*) 


+ a2(a&*)^ + a2(a5*)3 


where the four constants are fixed by the coordinates of the base 
points and the values of the curve slope in them. Thereby, 

3.5* 

Schlichting obtained for the Blasius profile higher values — — tl 

Ua 

those occiirring in figure 25 whereas the values of the slope are in 


than 


agreement. Actually, the curves 


-(a5*) in the center part of the 


(a6*) -region enclosed in each case will run somewhat higher or lower 
than indicated by figures 22 to 26. At any rate, however, they may be 
interpreted as a first approximation in the usual sense. 


These curves 


Pi 6* w 

— — over a5 represent sections through the 
^a 


"excitation mountain range" enclosed by the neutral curve as base curve, 
along the lines y^ = Const, or c^ = Const. By interpolation one 

obtains from them the maps of excitation represented in figiores 27 to 32 , 


in which the lines of equal excitation 


= Const*, can be interpreted 


as "contour lines'* of the "excitation mountain range," Instead of the 
lines of intersection c^ = Const., we plotted in figures 27 to 32 the 
lines 


UaRe^ 


= Const, 


the significance of which will be discussed later in section VII. Even 
at first consideration of .these excitation maps, a fundamental difference 
in the shape of the excitation mountain range is conspicuous according 
to whether the velocity profiles of the laminar friction layer lie in 
the region of decreasing pressure (p > 0) or of increasing pressure 

(3 < 0). 

In the region of decreasing pressure, the "excitation mountain 
range" has the form of a mountain with pronounced peak which is steeply 
ascending for a small Reynolds number Re*, slowly flattens after a 
larger Re* and shrinks to zero width and height for Re* — ><». The 
absolute height of the peak, that is, the maximum excitation increases 
with decreasing p. 


5A 


NACA TM 13i^3 


33 


In the region of increasing pressure, in contrast, the excitation 
mountain range changes behind the peak with growing Re* into a 
"mountain ridge" of constant width and constant contour profile. 

The properties of this "ridge" have been thoroughly Investigated 
by W. Tollmien (ref. 5). Since we made use of Tollmien’s theory, we 
must now briefly represent its results. 

Searching for a general instability criterion, W. Tollmien estab- 
lished that the frictionless perturbation equation (5) possesses for the 
laminar velocity profiles in the region of pressure increase - which 
have a point of inflection in contrast to those of the region of 
decreasing pressure - for Re*—*' » aside from the neutral solution 
existing for all profiles with the parameters 


a = 0 Cr = 0 cpnl = U (93) 

in addition, the neutral solution with the parameters 

a = ttg Cr = Ug cpnil = Ts (9^) 


the subscript s therein denotes the point of inflection. 

W. Tollmien calculated the eigenvalue ag and the eigenfunction 
cpg for the sinusoidal velocity distribution ( 30 ) which we also used 
from the frictionless perturbation equation (5) 


(aa) g = - p cot p 


<Ps = 


sin p I 


sin 


ps 


sin p 2^1 - .S-j 
(o < y < a) 


sin p aaf 1 
a 


(y>a) 


(95) 


( 96 ) 
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For this second neutral eigensolution, the phase velocity Cj. therefore 
equals the velocity of the basic flow at the point of inflection Ug. 

For the excitation and the circular frequency 3r in "the 

neighborhood of the neutral frequencies a6* = 0 and a&* = (a6*)g, 

W. Tollmien has derived the following formulas 


= n Hifll H^(a6*)3, Re* = 

5*2' Ua2 U„’6* 


where the subscript w signifies that the values have to be taken at 
the wall y = 0, and that 



with 


where 





(99) 


(100) 


We have calculated the "ridge contour profiles" of the excitation 
and the circular frequency according to the formulas (97) and (98) for 
the Hartree velocity distributions 3 = -0.10, -0.l4, -O.I6, -0,l8, 
-O.I98J the values for the derivatives of the velocity were not taken 
from the sinusoidal approximation distribution but from the exact 
velocity distribution. These "ridge contour profiles" of the excitation 
and the circular frequency have been plotted in figures 33 and 34. The 
curves corresponding to the formulas just mentioned have been drawn in 
solid lines) the transitions between the two curve arcs, interpolated 
somewhat arbitrarily, have been drawn in dashed lines. In the range 
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of very large Re*, the lines of constant excitation and constant cir- 
cular frequency were transferred from figures 33 and 3^ into figures 31 
and 32. 

Figure 35 represents the variation of the maximum excitation 


IhA 

\ Ua / max 

against the form parameter 3, separately for the excitation for finite 
Re* (peak of the excitation mountain range) and for Re* > 00 (ridge 
of the excitation mountain range). For 3 < -0.10 these two values 
seem more and more to approach one another so that with decreasing 3 
the peak becomes less and less pronounced, and the excitation starting 
from small Re* monotonically increases to the values of the "ridge 
contour profile." 

This figure and the preceding ones show clearly that the maximum 
excitation and accordingly the excitation in general is considerably 
larger in the region of pressure increase than in that of pressure 
decrease for smaller Re* as well. 

After thus having estimated the magnitude of excitation in the 
entire instability range of all velocity distributions occurring in the 
laminar friction layer, we shall discuss the physical conclusions 
resulting from our calculations for the position of the transition point. 


VII. DISCUSSION OF THE RESULTS 


Let the pressure distribution Ua(x) against the arc length of 
the cross section profile, and the oscillation of circular time frequency 
3j. be prescribed. Then this oscillation superposed on the boundary- 

layer flow travels downstream on a curve 


As we mentioned at the beginning, there pertains to every point x of 
the profile in the flow a fixed value of the Pohlhausen parameter \ 
which characterizes by way of approximation the velocity distribution 
in the boundary layer at this point. According to figure 3^ in refer- 
ence 3 one may coordinate to this parameter X. the Hartree parameter 3 
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and therewith one of the excitation maps calculated in the present 
report (figs. 27 to 32) . In order to obtain the excitation of the 
perturbation at a certain point x of the profile, one has there- 

36^ 

fore to read off the excitation (number of the "contour line") 

Ua 

on the corresponding excitation map at the point of the map determined 
by the pair of values of the Reynolds number 


Re* 


V 


g(x) 


and the dimensionless circular frequency 


U 


= f(x) 


a 


We shall start the discussion of the results of our above excitation 

Uoot 

calculations with the limiting case of a small Reynolds number Re = 

According to the explanations in section I (compare fig. 1) the stability 
limit lies, for small Reynolds numbers only, at the point where the 
laminar boundary layer separates. If the perturbation waves are so 
long that the curve 


B V 

* ■ 


intersects the instability region for the separation profile (3 = -O.I 98 ) 
it is very violently excited when entering this zone and leads quickly 
to transition to the turbulent flow pattern. The transition point then 
practically coincides with the separation point which we had denoted 
its rearward limit. If, on the other hand, the perturbation waves are 
very short so that the curve 


3pV 


f(x) 


does not intersect the instability region, the laminar layer separates 
without transition. 
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If the Reynolds number Re = increases, the stability limit 

shifts forward in the direction toweird the pressure minimum (fig. 1), 
the perturbation enters a region of instability further upstream for a 
larger value of the parameter 3 or the Pohlhausen parameter and 
is there initially excited to a degree which decreases the more the 
stability limit shifts forward but which increases due to the fact that 
the perturbation downstream (with decreasing 3 or X,) reaches insta- 
bility zones with rapidly increasing excitation. The transition point 
shifts frontward corresponding to the excitation which started earlier 
and is still strong. 

If we finally increase the Reynolds number Re = so that the 

stability limit shifts ahead of the pressure minimum, the excitation in 
turn starts accordingly soonerj however - and this must be regarded as 
the most important result of our calculations for the time being - in 
the region of decreasing pressure, the excitation is so slight that it 
generally attains amounts equalling the excitations produced in the 
cases treated just now only after having passed the pressure minimum. 

In this manner, one may easily give the theoretical explanation 
for the fact proved by many experiments, that the transition point even 
in case of very high Reynolds numbers rarely ever shifts ahead of the 
pressure minimum. Only in cases of very long acceleration sections 
and perturbations with very long waves where the perturbation does not 
too soon leave the instability region again, the limited excitation of 
the region of decreasing pressure will be sufficient to induce the 
transition still in the region of pressure decrease. However, these 
cases are rare in technical application. 

A detailed calculation of the degree of excitation which causes 
the transition is meaningful only in connection with corresponding 
experiments. Therefore, it will be postponed until these experiments, 
now in the preparatory stage, have been carried out. Probably one will 
have to regard the form of the pressure distribution in the region of 
pressure increase as the most important test condition since, according 
to our theoretical deliberations, the contributions to the excitation 
of perturbations in the region of pressure decrease are insignificant. 

The aim of this experimental investigation and of the excitation cal- 
culation to be performed simultaneously on the basis of the present 
report will be to find a connection* between the pressure-distribution 
form and the degree of excitation attained at the measured transition 
point so that it will be possible to calculate, inversely, the transition 
point for a prescribed pressure distribution from this relation. 
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VIII. SUMMARY 


As a contribution to the solution of the important, problem of the 
calculation of the transition point of a plane laminar flow, we had 
first determined (in an earlier report, according to Tollmien's method 
of small oscillations for the Hartree velocity distributions appearing 
in the boundary layer in case of decreasing and of increasing pressure) 
only the critical Reynolds number beyond which the perturbations super- 
posed on the laminar flow are excited. In connection with those cal- 
culations now, the excitation itself in the entire instability range 
of the perturbations was calculated. The excitation in the narrow 
instability range of decreasing pressure turns out to be very much 
smaller than the excitation in the more extensive instability range of 
increasing pressure; thus the known fact that the transition point 
generally does not shift ahead of the pressure minimum even in case of 
high Reynolds numbers may be explained on a theoretical basis, as shown 
in tables 1 , 2 , and 3 * 

Systematic experimental measurements of the transition point, 
together with calculations to be performed on the basis of the results 
given here, are to establish the connection between the variation of 
the pressure gradient and the degree of excitation which produces the 
transition and thereby a basis for determination of the transition point 
for prescribed pressure variation by calculation. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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TABUS 1 

SERIES CCBFPICIEWrS OF TEE SOIOTIOHS FOB THE BASIC FUW 
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5 COEFFICIEIITS OF THE 80UJTI0N8 <Pg* FOR THE BASIC FUW 
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^RIES COEFFICIEIITS OF THE SOLUTIONS FOR THE BASIC FLOW 

U * Up + (1 - Ufl) sin /3LI-® ABBREVIATION: p = ^ 

S' - 8 2 j 8lny2oCOBy2B 

Po* =0 Si* * - tan yggCl + ip) Sg* = ^ | tan?y2a(l + 2 ip) - 02^ 

P3* = - I tan y2g(l + Ip) - | tan3y2a(a^+ 3 ip) 

34* = ^ tan?y2a(l + 21p) + | tanV2s(l + 

P 5 * = - ^ tan yggCl + Ip) - ^ tankas (1 + 3ip) - ^ tan 5 y 2 s(l + 5ip) 

06* = tan^2s(^ + 21p) + ^ tanVaa^^ ^ 

®7 " ■ I^lo ^ tan^yggd + 3ip) - ^ tan^yggd + 5ip) - ^ tanTyggCl + Tip) 

3q* = ^ tanSyggd + 2ip) + SL tan^2Bd + ^tp) + ^ tan^2a(l + + ®iP) 


«o* “ ^ ®1* “ 5 ^ 20(1 + ^P) ® 2 * 

1 n^2 

©3* » - ^ tan y2a(l + ip) + ^ tan yggCl + ip) 

®‘‘* = lio - i t«n2yg,(l » aipl + gj 

<v,2nn 




6 6 


*■* = i t«n-=yg,(l » aipj + ^ 

*5^ * ^ ^28^^ * - I ♦ 31p)j^ og*lt«ii yggd - Ip) 

*6* = - 5^ ♦ i t“^a8(i ♦ 2ip) - ^ »pVa.a ♦ *‘ip)] - j^ji + ^ t«t^S8(i ♦ 2ip^ ♦ 

»_* s _ __JL 'f^n V— ("X A, -f-n^ X **2. .r^T v— /1 ^ 4'n^ ^ 29 •*>.»^3u' ^ ^ 1 /l . ^j[p^ . 


■ 5 ^ "■ I680E ^ _, ioooL- =(U 

tan yggd ♦ ip) + tan y 2 a<l ♦ ^P) * ^ tanSy^gd + 31 p) + ^ tanSy^gd + 

4 r ^ 6 

tan yggd * Ip) - tanSygad t 31 p^ t tan yggd + Ip) 

fsLli * 4 . 4 . §57 4,„1l,- m J. Jl'fTk\ 4. _3 -fan 


b2* 

b,* 


J 2 bi»ei 00 

*8* ' 3 ^i §0 * !^|| * ^ tai^yagd + aip) ♦ ^ tanV 2 e(l + *‘lP) ♦ ^ tan8ygg(l + 6lp^ + 

^|T - ^ tan^yggd * aip) - tanVagd * 4lp^ - ^|l t Ji tan?yggd . 2 ip^ . ^ 

- H + tan^yggd + aip^ + | <«a® 

tan yag(l ♦ Ip) - | tan 3 ygg(l + 31^ - ^ tan ygg(l + Ip) 

- ^ * 5 tnn^ggd + 21p) + i tanVaad + - “a^p + ^ tan2ygg(l + aip^ + ^ 

** 5 * ' ■ [iTO t ^P) t tanSggd + 51 py + 02®j^ tan ygg(l + Ip) - tan 3 ygg(l + 3 lp^ . 

Og^ 

^ tan yggd + ip) 

' 1} ^ tai^ggd t aip) + ^ tanVa,(l ♦ **lp) + ^ tan^yggd + 6lp^ + 

o»2[X ♦ . m x. tan2v„.ri + ainTl + ^ 




__ ' — 

- tan yggd + Ip) - ^ tan3y2ad ^ 3ip) - ^ tan^y^^d + 5ip) - ^ tan'^ygsd + Tlp) 


' iSiTO 


fV/ — 

_- a^ioo yp8<^ * * 33^^00 **“^28 <i * 3 ^P) ♦ i^ooo * 5 ipj t 

: 1 ^ tan yggd t Ip) * tan 3 ygg(l ♦ 31 p^ - ^ tan yggd . Ip) 

^ tan2yggd ♦ 21p) + ^ tan^BaCl + tip) + ^ tan^yggd + 6lp) + 

5^ tanSyggd ♦ 8lp) + og2 Q * 21P) - tanVagd ♦ Hp) - 

ieiwioo t 6lp^ t ogt|^ ♦ tanSyggd . aip) - ^^^5 *“'^28(1 * ‘‘^P^ * 

“a® Q iok - wife ^ 
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TABLE 2 


SERIES COEFFICIENTS OF TEE SOLUTIONS , qpg* FOR THE BASIC FLOW 


U = Ug (1 - Ue) sin 5) ABBREVIATION: p = 

\S, - B ei 0ln 


yga COB ygg 


- Concluded 


^28 82 = + I ^ 83 = “ I ^28 ^ I ^28 


+.«n2^ 

24 




4 


lio =^2s * 


fo( - ="2 b - ^ ^28 

* ^ tanSyg, * - f tan^ya, + | tanSr^^j . . 11 . M tan2y, 

2 / X 

■ 5^ ’'SB + ^p8 ’'^S - ^ tan 3 ygj - | tan^yagj + 

(- ^ ’'28 - tan 3 ya,J . ^ 0^6 tan y^, 

• lok ^“^2" * Wo\^ ^ * is *‘"^’'2» * is **"‘^28 ♦ II tanfiya^ t 
sisoj^ - ^ '^‘"®’'2 b ♦ 7m *““*^2^ * 15 ^ [[_®^ * ^ tanSyg^ t ^ 


oo6 

28 I * Tio 


0 h^ = - 1 


- 2 tan^ygg + 02 ^ ^ = - J ^28 " ^ tan^ygg - i og^ tan ygg 


i + A tan2yg^ - i tanVag ’ ^ “a®^^ t f tan2yg,j t i 


1 4 


59 ^ 1 

tan yoo + — 
288 2 b 3 


L da** tan yag 


U " 25 , ^ 

+ i tanBygg + ^ tanSyag + ^ oa® ^ - 227 tan yag + tanSyg^j - i og' 

■ iie • '^"*’'28 - ^ t«Va8 - ^ tanfiyag + ofiSj^l + ^ tan2ygg - ^ tanVa^- 
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- • 5^ ^ “" ='2* 2 ^ 
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25920 ’'28 
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324000 7 
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TABLE 3 

VALUES OF THE DIFFERENTIAL QUOTIENTS OF THE REAL AND OF THE IMAGINARY 
PART OF THE FUNCTION WITH RESPECT TO ' Tio 


RO 

8 f^ 

^1^0 

8 Fj_ 

^Ro 

^0 

8 f^ 

8 Fj_ 

^RO 

-2 

0.135 

-0.226 

1 

-3.6 

0.396 

-0.066 

-2.1 

.121 

-.230 

-3.7 

.399 

-.030 

-2.2 

.117 

-.235 

- 3.8 

.395 

-.010 

- 2.3 

.118 

-.239 

-3.9 

.379 

-.055 

- 2 .k 

.121 

-'.242 

- 4.0 

.360 

-.100 

- 2.5 

.128 

-.245 

■ - 4.1 

.330 

-.155 

- 2.6 

.135 

-.246 ■ 

- 4.2 

.297 

-.205 

- 2.7 

.145 

-.246 

-4.3 

.250 

-.270 

- 2.8 

.157 

-.243 

- 4.4 

.180 

-.338 

- 2.9 

.172 

-. 237 ' 

- 4.5 

.080 

-.341 

- 3.0 

.191 

-.226 

- 4.6 

.020 

-.282 

- 3.1 

.219 

-.210 

- 4.7 

-.026 

-.226 

- 3.2 

.256 

-.188 

- 4.8 

-.045 

-.182 

- 3.3 

.296 

-.162 

-4.9 

-.057 

-.150 

- 3.4 

.340 

-.133 

-5.0 

- . 062 

-.118 

- 3.5 

.380 

-.102 

1 

1 
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Figure 1 .- Frontward limit (stability limit) and rearward limit (separation 
point) for the transition point on a cylindrical body. 
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Figure 13 


Figures 11, 12, and 13.- Polar diagrams for determination of the curves of 

constant excitation. 
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Figures 17 and 18.- Polar diagrams for determination of the curves of con- 
stant excitation. 
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Figure 24 

Figures 22 , 23, and 24.- The excitation — as a function of the reciprocal 

2TT 

perturbation wave length ~ ~ ^®r constant critical velocity c^. 
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Figure 25 
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Figures 25 and 26.- The excitation 

27T 

perturbation wave length ^ ^ ~ 


— — as a function of the reciprocal 
for constant critical velocity c^. 
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few laminar boundary-layer profiles. 




few laminar boundafy-layer profiles. 
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Figure 31.- The curves of constant excitation in the instability regions of a 
few laminar boundary-layer profiles. 
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Figure 35. 





/Pi5*\ 

- The maximum excitation ) as a fxinction of the form 

\ *^a /max 

parameter p of the laminar boimdary -layer profiles. 
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